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We present here three analyses of semileptonic Y>-meson decays based 
on the 2.92 fb“^ of data collected by the BESIII experiment in 2010 and 
2011 at the p(3770) peak. For the decay —?• A^vr+e+r'e) hs branching 

fraction is measured over the whole rriK-K region and in the K (892)° win¬ 
dow, respectively. A partial wave analysis (PWA) is performed, indicating 
an (S-wave contribntion of about 6%. The S'-wave phase and the form fac¬ 
tors are measured by the PWA and in a model-independent way. For 
the decay —?• uje'^Ve, an improved measurement of the branching frac¬ 

tion is performed and the form factors are determined for the first time. 
D+ —>■ (pe'^i'e is searched and an improved upper limit at 90% confidence 
level is set. For the decay —>■ K^e'^ne, its branching fraction is mea¬ 

sured for the first time and the CP asymmetry is presented. The product 
of the hadronic form factor and the CKM matrix element, ff^{0)\Vcs\, is 
also determined in this decay. 
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1 Introduction 


The analyses presented in this paper are based on the 2.92 fb“^ of data sample [T] 
collected in annihilations at the 'i/^{3770) peak, which is accumulated with the 
BESIII detector [2] operated at the double-ring BEPCII collider. At the 'ijj{3770) 
peak, D mesons are produced in pairs, thus allowing us to use the technique of 
tagging D-meson decays |3]. If a decay of one D meson (“tagged decay”) has been 
fully reconstructed in an event, the existence of another D decay (“signal decay”) 
in the same event is guaranteed. In these analyses D~ is hrstly reconstructed in 
one of the six decay channels: D~ —> A'+7r“7r“, D~ A'+7r“7r“7r°, D~ 

D~ Kg7r~n^, D~ —>■ A'^7r“7r“7r+, and D~ , then semileptonic signals 

are searched at the recoiling side. Charge conjugated decays are implied similarly. 

2 Study of K~77^e^i'e 

The invariant matrix element of the decay of —?• is the product 

of a leptonic and a hadronic current. Thus it’s convenient for us to study the Kn 
system in the absence of interactions with other hadrons and the hadronic transition 
form factors. In this analysis we measure the different Kt^ resonant and non-resonant 
amplitudes that contribute to this decay, which can provide helpful information for 
the 5-meson semileptonic decays. Also we measure the dependent transition form 
factors, where is the invariant mass square of the lepton pair of the eVg system. This 
can be compared with hadronic model expectations and lattice QCD computations |1]. 

The signal yield is obtained by simply counting since the background level is lower 
than 1%. Using the tagging technique introduced in Sec.[^ the branching fractions are 
measured to be (3.71±0.03±0.08)% over the full rrixTr range and (3.33±0.03±0.07)% 
in the K (892)° region, respectively. 

The 4-body De 4 decay can be uniquely described by five kinematic variables |3]: 
Kn mass square (m^), ei^g mass square (q^), the angle between the tt and the D 
direction in the Kt: rest frame {Ok), the angle between the Ug and the D direction 
in the ei^g rest frame (6*e), and the angle between the two decay planes (y). Its PDE 
can be expressed based on the five variables and then the PWA is performed using 
an unbinned maximum likelihood method. 

The PWA shows that besides the dominant K (892)°, an ^-wave contribution 
making of the non-resonant background term and the 770(1430)° accounts for (6.05 ± 
0.22 ± 0.18)%. Other components can be neglected. The K (892)° parameters 
are determined: = (894.60 ± 0.25 ± 0.08) MeV/c^, ~ (46.42 ± 

0.56 ± 0.15) MeV/c^, and the Blatt-Weisskopf parameter tbw = 3.07 ± 0.26 ± 
0.11 (GeV/c)“h The parameters defining the hadronic form factors are also mea¬ 
sured: rv = ^ = 1.411 ± 0.058 ± 0.007, ra = ^ = 0.788 ± 0.042 ± 0.008, 
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mv = (1.8lla?7 ± 0.02) MeV/c^, ttia = (2.6llg;?^ ± 0.03) MeV/c^, ^i(O) = 0.585 ± 
0.011 ± 0.017. my is firstly measured for this decay. Corresponding projections over 
the five kinematic variables are illustrated in Figure 
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Figure 1: Projections onto each of the kinematic variables, comparing data (dots with 
error bars) and signal MC weighted by PWA solution (line), assuming that the signal 
is composed of the ^-wave and K (892)°. The shadowed histogram is the estimated 
background. 


In the above PWA process Ss depends on rriK-K according to the LASS parame¬ 
terization. Then we measure 5s in a model-independent way. We divide the rriK-w 
spectrum into twelve bins and perform the PWA fit with 5s in each bin as twelve 
additional fit parameters (within each bin the phase is assumed to be constant). Fig¬ 
ure illustrates the comparison of the model-independent measurement with that 
based on the LASS parameterization. 

In the PWA, the helicity basis form factors are assumed to depend on according 
to the spectroscopic pole dominance (SPD) model. In order to achieve a better under¬ 
standing of the semileptonic decay dynamics, we also measure the form factors in a 
model-independent way using the projective weighting technique, which is introduced 
in Ref. [6]. 

Only candidates in the W (892)° region ([0.8,1] GeV/c^) are used, so that we 
can neglect other components than the K (892)° and the non-resonant A-wave. The 
decay intensity of D~^ —?• K~Tr~^e~^Ue can be parameterized by form factors describing 
the decay into the vector meson K (892)°: H_{q‘^,m), Ho(q^,m), and by 

an additional form factor ho{q‘^,m) describing the non-resonant S'-wave contribution. 
The form factors are measured by weighting the distributions based on the angular. 
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Figure 2: The 5'-wave phase variation versus ttikw assuming that the signal is com¬ 
posed of the ^-wave and K (892)°. Model-independent measurement (points with 
error bars) is compared with the result based on the LASS parameterization (solid 
line, la deviation is marked by dashed line). 


The results are shown in Figure They are consistent with the SPD model with the 
parameters obtained from the PWA, and with the results reported by CLEO-c [7j. 





Figure 3: Form factors measured in this work (squares) compared with the CLEO- 
c results (circles) and with the PWA solution (curves). Error bars represent both 
statistical and systematic uncertainties. 


3 Measurement of the form factors in the decay 
D+ ^ cje+z/g Sind search for the decay 

The decay uje^v^ has similar dynamics as K *Neglecting 

the mass of the electron, the transition matrix element of D to vector meson can 
be decomposed into contributions from one vector V{q^) and two axial-vector (Ai, 
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A 2 ){q‘^) form factors, where is the invariant mass square of the system. A 
precise measurement of the branching ratio and the form factors provide opportunities 
to test the standard model and the theoretical calculations. 

The decay —?• has not been observed at present. The (p (ss) has different 

quark composition from the D meson (cd), so the process can only proceed either 
through u} — (f) mixing or non-perturbative “weak annihilation” (WA). A measurement 
of the branching fraction can discriminate which process is dominant. 

The signal yields are determined by the variable U, the difference between the 
missing energy and momentum jS]. The yield of D~^ uje'^Ue is obtained from a fit 
to the U distribution as shown in the left plot of Figure]^ And that of —>• 0 e+z/e 
is obtained by counting the number in the signal region [-0.05, 0.07] GeV as shown 
in the right plot, indicating no significant excess of signal events. The results are 
concluded in Table [T| which are improved compared with previous reports. 




Figure 4: Left: fit (solid line) to the U distribution in data (points with error bars) 
for D~^ —>■ cue'^i'e- The total background is shown by the filled curve, with the peaking 
component shown by the cross-hatched curve. Right: the U distribution for D~^ 

in data (points with error bars) and signal MC with arbitrary normalization 
(solid histograms). The arrows show the signal region. 


Table 1: Measured branching fractions and a comparison with the previous measure¬ 
ments. For the first uncertainty is statistical and the second systematic. 


Mode 

This work 

Previous [8l [9] 


(1.63 ±0.11 ±0.08) X 10“^ 
1.3 X 10-5 (90%C.L.) 

(1.82 ±0.18 ±0.07) X 10-5 
9.0 X 10-5 (90%C.L.) 


In order to measure the form factors in the decay a five-dimensional 

maximum likelihood fit is performed in the space of cos 6 *i, cos 6*2 and y, 

whose definitions are similar with those described for the decay —>■ A'7re'''i^e in 
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Secj^ The form factor ratios are determined from the fit: ry = 1-24 ± 0.09 ± 0.06, 
r 2 = 1.06 ±0.15 ±0.05, which are measured for the first time in this decay. The fitted 
projections over the five variables are illustrated in Figure 



Figure 5: Projections of the data set (points with error bars), the fit results (solid 
histograms) and the sum of the background distributions (filled histogram curves) 
onto (a) rri^, (b) (c) cosdi, (d) cosd 2 and (e) y. 


4 Study of decay dynamics and CP asymmetry in 

decay 

In charged D meson decays, CP asymmetry occurs when the absolute value of the 
decay amplitude for D decaying to a final state is different from the one for the 
corresponding CP-conjugated amplitude. The most optimistic model-dependent es¬ 
timates put the SM predictions for the asymmetry as 0(10“^) or below (TU], so an 
observation of any CP-violating signal would be a sign of new physics. The hadronic 
matrix element of D to pseudoscalar meson process can be decomposed into contri¬ 
butions from longitudinal and transverse form factors. Neglecting the lepton mass, 
only the transverse form factor contributes. Various theoretical techniques provide 
slightly different dependencies of the form factors, High precision measurements 
of the partial decay width over different ranges of (f‘ will distinguish which method 
correctly describes the non-perturbative dynamics of QCD. The —>■ decay 

is investigated in this work with its branching fraction and CP violation firstly mea¬ 
sured. dependence of the form factors are measured based on different theoretical 
models for the first time as well. 
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In the process of searching for D+ —> K^e^Vf, signals, the basic idea for 
search is explicitly introduced. Most will penetrate the MDC, interact with the 
scintillators in the EMC and provide the position information. With all charged 
particles in the event reconstructed and applying energy-momentum conservation to 
all the particles, the direction of the missing can be determined. We then examine 
whether there’s a good neutral shower deposited in the EMC along the predicted 
direction. 

The signal yield is obtained by fitting the distributions of the beam constrained 
mass (Mbc) of the corresponding tagged D candidates, as illustrated in Figure 
Branching fractions of the six tag modes are measured separately for D~^ and D~. 
Then the weighted branching fractions for D+ and D~ are obtained: B{D~^ —)■ 
Kle+Ue) = (4.454 ±0.038 ±0.102)%, B{D- K^e-u^) = (4.507 ± 0.038 ± 0.104)%. 
The combined averaged value is 

B{D+ Kle+Ue) = (4.481 ± 0.027 ± 0.103)%, (1) 


which agrees well with the measurement of B{D^ —?• Kge~^Pe) by CLEO-c [H]. The 
CP asymmetry is determined to be 


_ B{D+ ^ K^e+ue) - B{D- ^ i^^e-^) 
“ B{D+ Kle+Ve) ± B{D- K^e-Ve) 


(-0.59 ± 0.60 ± 1.48)%, 


( 2 ) 


which is consistent with the theoretical prediction in Ref. ca 



Figure 6: Fits to the Mbc distributions of the tagged D~ candidates for data, requir¬ 
ing —7- Kge~^Pe. The dots with error bars are for data, and the blue solid curves 

are the results of the fits. The green dashed curves are the fitted backgrounds. 
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We perform simultaneous fits to the distributions of the observed e+z/e 

candidates as a function of to determined the product of the hadronic form factor 
and the CKM matrix element f^{0)\Vcs\- The form-factor shape is described based on 
the simple pole model (mpoie), the modified pole model («), two-parameter series ex¬ 
pansion (ri), and three-parameter series expansion (ri, r 2 ). As an example, Figurej^ 
shows the simultaneous fits using the two-parameter series expansion model, corre¬ 
sponding to the results: f^{0)\Vcs\ = 0.728±0.006±0.011, ri = -1.91 ±0.33±0.24. 



Figure 7: (Color online) Simultaneous fits (blue solid curves) to the numbers of 
data (points with error bars) as a function of with the two-parameter 
series expansion parameterization. The red dashed curves show the signal, while the 
violet, yellow, green and black curves refer to different kinds of backgrounds. 


5 Summary 

Based on 2.92 fb“^ of data collected by the BESIII experiment three semileptonic 
decays are analysed. In the study of D+ —)■ A'“7r+e+zze, its branching fractions are 
measured: B{D^ K~7T~^e~^U(,) = (3.71±0.03±0.08)%, B{D^ K~Ve)[o.%,i] = 
(3.33 ± 0.03 ± 0.07)%. A PWA is performed and an ^-wave contribution is found to 
account for (6.05±0.22±0.18)%. The S-wme phase and the form factors are measured 
both by the PWA and in a model-independent way, showing good consistency. For 
the decay —>■ the branching fraction is measured with a higher precision: 

B{D^ uje^Ve) = (1.63±0.11 ±0.08)%. Its form factors are determined for the first 
time: ry = 1.24 ±0.09 ±0.06, r 2 = 1.06 ±0.15 ±0.05. The rare decay H"'' —?• is 
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searched for and an upper limit at 90% confidence level is set: ^e+r'e) < 1-3x 

10“^. For the decay —?• Ki^e'^Ug, we present the first measurement of the absolute 

branching fraction ^ K^e'^p^) = (4.481±0.027±0.103)%. The CP asymmetry 
is also determined Acp = (—0.59 ± 0.60 ± 1.48)%, indicating no significant CP 
violation. The product f^{0)\Vcs\ is measured based on three theoretical form-factor 
models. Taking the two-parameter series expansion parameterization, f^{0)\Vcs\ = 
0.748 ±0.007 ±0.012. 
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